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Neutrino-neutrino interactions in dense neutrino streams, like those emitted by a core-collapse su- 
pernova, can lead to self-induced neutrino flavor conversions. While this is a nonlinear phenomenon, 
the onset of these conversions can be examined through a standard stability analysis of the linearized 
equations of motion. The problem is reduced to a linear eigenvalue equation that involves the neu- 
trino density, energy spectrum, angular distribution, and matter density. In the single-angle case, 
we reproduce previous results and use them to identify two generic instabilities: The system is sta- 
ble above a cutoff density ("cutoff mode"), or can approach an asymptotic instability for increasing 
density ("saturation mode"). We analyze multi-angle effects on these generic types of instabilities 
and find that even the saturation mode is suppressed at large densities. For both types of modes, 
a given multi-angle spectrum typically is unstable when the neutrino and electron densities are 
comparable, but stable when the neutrino density is much smaller or much larger than the electron 
density. The role of an instability in the SN context depends on the available growth time and on 
the range of affected modes. At large matter density, most modes are off-resonance even when the 
system is unstable. 

PACS numbers: 14.60. Pq, 97.60.Bw 
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I. INTRODUCTION 



Neutrino flavor oscillations in a supernova (SN) are 
strongly suppressed by matter effects |l| until the neutri- 
nos pass through the usual MSW region 2-5] far out in 
the envelope of the collapsing star. However, neutrino- 
neutrino interactions [y, 0] , through a flavor off-diagonal 
refractive index, can trigger self-induced flavor conver- 
sions [a-tl3| . This collective effect tends to occur between 
the neutrino sphere and the MSW region and can lead 
to strongly modified neutrino spectra, showing features 
such as spectral swaps and splits [lJ-Qjj ; for a review see 
Ref. [20] . The overall scenario, supported by heuristic ar- 
guments and numerical examples, is that deep inside the 
SN core, the system performs "synchronized oscillations" 
with an extremely small amplitude, i.e. every neutrino 
remains essentially stuck in its initial flavor cigenstatc. 
As the neutrinos stream outwards, there is a sharp onset 
radius where "bimodal" oscillations begin: Some ranges 
of modes start pendulum-like oscillations [ll|, l2lM23| , ex- 
changing their flavor content with each other without af- 
fecting the flavor content of the overall system. 

This scenario engenders a crucial simplification for the 
treatment of neutrino transport in SN simulations. At 
high densities, where neutrinos collide frequently, it is 
enough to solve the transport equations for each flavor 
separately, ignoring oscillations entirely. On the other 
hand, flavor conversions at larger distances can be treated 
ignoring neutrino collisions and absorption, i.e. as a pure 
propagation problem. So collisions and flavor oscillations 
are phenomena that are assumed to be taking place in 
different regions of the star and can be treated indepen- 
dently. When the radial distance where bimodal oscilla- 
tions begin is far away from the SN core, this assumption 
is valid and the flavor conversions do not affect the SN 



dynamics. Recent studies dedicated to the SN accretion 
phase, under simplifying assumptions, once more confirm 
this picture (2^.|25J. 

However, what is missing is a systematic approach to 
decide, without solving the equations of motion, if self- 
induced flavor conversions occur for given neutrino spec- 
tra (flavor-dependent energy and angular distribution), 
overall neutrino density, and matter density. Formal sta- 
bility criteria exist only in the "single-angle approxima- 
tion" where it is assumed that all neutrinos feel the same 
neutrino- neutrino refractive effect. In this case the ana- 
lytic pendulum solution has been found and its existence 
and parameters can be calculated from the neutrino spec- 
trum and density alone [19(. 

On the other hand, the current-current nature of the 
low-energy weak-interaction Hamiltonian implies that 
neutrinos in the background of an anisotropic neutrino 
flux experience a refractive effect that strongly depends 
on direction. For some energy spectra, these "multi-angle 
effects" have little impact, whereas in other cases they 
completely change the solution. A SN neutrino flux with 
a vanishing net v e flux is unstable everywhere and shows 
quick multi-angle decoherence (26[ ■ If the v e flux is large 
enough, this effect is avoided if the energy spectrum is 
simple [27|. More complicated spectra can be unstable 
even for a large v e flux at any density in the single-angle 
case, but the instability is suppressed by multi-angle ef- 
fects [28|, (29[ . In addition, the presence of ordinary mat- 
ter causes a multi-angle suppression of the bimodal in- 
stability [301 . On the other hand, it was claimed that 
for nontrivial angular distributions, as may exist when 
different flavors are emitted from very different neutrino 
spheres, there is a novel multi-angle instability [12] ■ De- 
viations from the usually assumed cylindrical symmetry 
may also have an important influence [13j. 



Although our problem is nonlinear and therefore would 
seem intractable, noting that an instability must occur 
in order for the onset to take place leads to a surprising 
simplification. In the dense SN matter well inside of the 
MSW region, the matter effect is so large that neutrino 
propagation eigenstates are essentially identical with fla- 
vor eigenstates. This means that in the weak-interaction 
basis, the flavor matrices of occupation numbers are al- 
most perfectly diagonal, allowing us to linearize the equa- 
tions of motion (EoMs) in terms of the small off-diagonal 
elements. An instability is equivalent to some of these 
small elements starting to grow exponentially. All we 
need to do is linearize the EoMs, perform a Fourier anal- 
ysis, and seek exponentially growing solutions of the rel- 
evant eigenvalue equation. This is simple even in the 
multi-angle situation. Almost certainly this approach 
can be extended to cases without the usually assumed 
cylindrical symmetry around the radial direction. 

Studying the stability of a strongly coupled and non- 
linear system in the small-amplitude limit is a standard 
technique. In the present context it was put forth in 
Ref. [l3|. However, the method was carried only to 
schematic cases of a small number of neutrino momen- 
tum modes, leaving open how to apply it to realistic sit- 
uations. 

While a stability analysis provides crucial insight, we 
stress that it alone is not enough to assess the impact in 
the SN context. Only a small range of modes may partic- 
ipate in the bimodal oscillation or the growth rate may 
be too small on the available time scale. (Of course, both 
the growth rate and the location of the resonance on the 
energy-angle spectrum are found by solving the eigen- 
value equation.) For example, the classic single-angle 
case with an initial Fermi-Dirac spectrum of only v e and 
v e is always unstable and the usual concept of a syn- 
chronization radius does not apply. However, for large 
neutrino densities, the growth rate is small and only a 
narrow range of infrared modes is affected. A "visible" 
effect arises only at a quasi-onset radius where the growth 
rate begins to compete with the overall evolution time 
scale and the resonance begins to move into the main 
part of the spectrum. We will here largely avoid such 
issues and concentrate on setting up the method and dis- 
cussing simple but informative schematic cases. Apply- 
ing this method in a realistic SN context will be left to 
future work. 

Our study in Sec. |TT] leads to the linearized equations 
of motion at a large distance from the neutrino source 
in the two-flavor case, with an azimuthally symmetric 
neutrino emission. In Sec. IIII1 we present the stability 
analysis in the single-angle approximation, and illustrate 
it in Sec. IIVI with the examples of box spectra where 
the results may be understood analytically. In Sec. |V] 
we point out some special features of realistic spectra 
that do not vanish at low energies. Sections IVII and IVIII 
demonstrate how the single-angle results are modified by 
the inclusion of multi-angle and matter effects. The lat- 
ter one also analyzes a realistic SN spectrum using the 



insights obtained from the box spectra. In Sec. IVIII1 we 
show that a multi-angle spectrum with a zero crossing is 
unstable in both hierarchies if the lepton asymmetry is 
small. In Sec. IIXI we conclude with a brief summary of 
our findings and an outlook on future directions. 



II. EQUATIONS OF MOTION 

A. Effective Hamiltonian 

We write the EoMs in terms of n-p x rip matrices of 
occupation numbers [TJ, l3l| , where n? is the number of 
flavors. We denote these matrices by Qe,v, where the 
velocity vector v with |v| = 1 describes the direction 
of motion and the energy E is taken to be positive for 
neutrinos and negative for antineutrinos. Also, while the 
diagonal entries are equal to the occupation numbers for 
neutrinos, they are the negative occupation numbers for 
antineutrinos. In the context of a Boltzmann collision 
equation for mixed neutrinos, one uses positive occupa- 
tion numbers in both cases and describe each mode by its 
momentum p Q. Our choice of signs, however, allows us 
to include neutrinos and antineutrinos on the same foot- 
ing and we will never have to distinguish between them: 
The antincutrino spectrum is simply a continuation of the 
neutrino spectrum to negative energies. In the language 
of flavor polarization vectors, our convention agrees with 
the neutrino flavor isospin construction |14l |. 

The EoMs for the time evolution in a homogeneous 
medium arc 



idtQE.v — [H£,v,te,v] 



(1) 



We use sans-serif letters for matrices in flavor space. The 
Hamiltonian matrix is 



M 2 r- 
Hb.v = — + V2G F Ni 



(2) 



\/2G F 



E' 2 

dv'-— ^ v ,(I-vv') 



where M 2 is the neutrino mass-squared matrix and ISh 
the matrix of net charged-lepton densities which in the 
flavor basis is ISh = diag(n e — n Sl n^ 



np,,n T ~nr , 



In an 



isotropic medium, the v • v' term drops out and the 
neutrino-neutrino term has the same structure as the 
matter term: The phase-space integral over qe, v amounts 
to the difference between neutrino and antineutrino den- 
sities. In the presence of macroscopic matter fluxes, 
there would also be current contributions in the matter 
term 13011 . 



B. Azimuthal symmetry 

Henceforth we assume azimuthal symmetry around 
some preferred direction, usually the radial direction in 



the SN case. The azimuthal integration provides 
1 — v • v — > 1 — v ■ v , 



(3) 



where v and v' are the components of v and v' , respec- 
tively, along the symmetry direction. Thus v = cos i? 
with $ the trajectory angle relative to the symmetry di- 
rection. 

Following the Appendix of Ref. j27| we consider the 
EoMs expressed in terms of the radial coordinate. We 
introduce an arbitrary sphere with radius R that we call 
neutrino sphere where we specify the inner boundary con- 
dition for neutrinos that are assumed to stream only out- 
ward. Every angular mode is described by its emission 
angle $r relative to the radial direction at that sphere 
(Fig. 3 of Ref. 20]) in terms of the variable u = sin 2 ## 
which lies in the range < u < 1. The u variable has 
the property that the modes are uniformly distributed 
on < u < 1 if the emission at the neutrino sphere is 
isotropic into space in analogy to blackbody emission. 

At radius r, the radial velocity of a mode with angular 
label u is 



R 2 



V u ,r = \/l-^U. 



In analogy to Ref. |27( we introduce the matrices 

r 2 E 2 



<t> 



E.u.r 



2n 



QE,u,r 3 



(4) 



(5) 



where we have included a factor Airr 2 , so that the inte- 
grated quantity 



<t>, 



dE I du <$>E.u r 



(G) 



represents the flux through a sphere of radius r whose 
trace is conserved. 

The EoMs for the flux matrices as a function of radial 
coordinate are 



id r Q>E,u,r — [H£, Uir , <£>£,„,,.] 
with the Hamiltonian 

V2G F N e ) — 



(7) 



H, 



2E 



*7T' Jo \ v u,r v u',r J 

where u , r = f_™ dE<& E .u,r- 

C. At a large distance from source 

We are interested in the evolution far away from the 
neutrino sphere where the flavor conversions are expected 
to begin. Therefore, we use the expansion 



«u,r = l 



u R 2 



(9) 



Moreover, we introduce the dimcnsionless matrices L = 
Ni/(n e — rig) and FE.u,r = ^e.u^/^pA^-)- Note that 
we normalize the charged-lepton density to the local net 
electron density, whereas the neutrino flux matrices are 
normalized to the total v e flux $ Pe (R) at the neutrino 
sphere. If we use the flavor basis, with these normaliza- 
tions we have L ee = 1 and J_ dE J Q du F^ur = ~ 1 f° r 
all r where oscillations have not yet begun. 

We also introduce the coefficients with dimension of 
inverse energy 



A r = V2 G F [n e (r) - n s (r)} 

\/2GF$ Pe (i?) 
^ R ~ 4^R 2 ' 

In terms of these coefficients, we have 



id r Fi 



[H 



E,u,r i *u,r\ 3 



(10) 



(ii) 



with 



lE.u.r 



M 2 
2~E 



= ^ + A r L 1 + ;^ 



uR 2 
2 7 2 



Mfl- 



R 4 



du' ■ 



:i2) 



as the Hamiltonian at the lowest-order in (R/r), with 
F ur = f dEFE,u,r- The first line on the right hand 



side of Eq. (|T2j) is the "vacuum plus matter" Hamiltonian 
Hj3 «V at wnne the second line is the neutrino- neutrino 
Hamiltonian H^" . 



D. Two-flavor case 

For the rest of this paper, we restrict ourselves to the 
two-flavor scenario, with flavors e and x, and we intro- 
duce the variable ui = \Am 2 \/2E 7 the vacuum oscillation 
frequency, to describe the different modes. In the con- 
text of flavor oscillation physics, w is a much more natural 
variable to describe the neutrino spectrum than the en- 
ergy E. Note that since E is taken to be negative for 
antincutrinos, they are represented by negative u> values. 
Since the trace of the Hamiltonian does not contribute 
to the time evolution, we write 



M 2 

2E 



co_ / cos26» sin26> 
2 l-sin20 -cos 26* 



A r L = -r- 



1 
-1 



(13) 



in the flavor basis, after removing a term proportional 
to the unit matrix. We take the mixing angle to lie 
in the first octant < 9 < ir/4. In this case the +(— ) 
sign stands for inverted (normal) mass hierarchy. In the 
following discussion, we shall consider inverted hierarchy. 
For obtaining results with normal hierarchy, we will have 
to multiply the u> term by a factor of —1. 



The flux matrices at the neutrino sphere are 



LU,U,R 



UJ,U 





(14) 



where the (j)ui,u are differential fluxes in the variables w 
and u. The normalization of F used here implies that 
J_ du> L duifiu u — —1. Note that (f> u ,u for antincutri- 
nos (lu < 0) corresponds to the negative of their occupa- 
tion numbers. Finally, in the flavor basis we write 



TrF, 



S 



(15) 



where g^^ = 4>tj u ~ 4>Z u i s the usual difference spec- 
trum, except that it is now also differential with regard 
to the direction variable u. The initial conditions at the 
neutrinosphere for the Hermitian matrix S WiMj7 . are 



~^UJ,U. 



1 
-1 



(16) 



It satisfies the EoMs 

with the neutrino-neutrino part of the Hamiltonian 



H ^u.r = Mr / du' (U + U') 



and 



+ OC 



2r 4 



dw'^-S w>v . (18) 



(19) 



The effective neutrino-neutrino interaction energy de- 
clines with r -4 . 



Next we write the S matrices in components in the 
flavor basis 



S 



a* _„ 



(22) 



where s^^^ is the r-dependent swap factor. It specifies 
how much the flavor content of the given mode has been 
swapped relative to the initial condition. We have the 

= 1. Likewise, 



normalization s^, ur + \S u , u ,r 



H 



UJ,U,r — I TT* 7 



(23) 



Then the EoM for the off-diagonal component is 

The components of the Hamiltonian matrix are explicitly 

LU + U X r 



*luj'U,r 



H... 



-^- du'(u + u')l du' g u 'u' Su^u', 
-^- du'(u + u')l du' g u > u > S u >, u > 

Z JO J-oo 



(25) 

In the absence of all interactions, the rotation-averaged 
EoM is 

implying the free precession solution 

c* - -iu(r-R) Q „ (07\ 

F. Small-amplitude expansion 



E. In a co-rotating frame 

We go to a rotating frame where the common mat- 
ter term drops out and where the vacuum term oscil- 
lates quickly, averaging the off-diagonal term to zero [3 Of . 
Moreover, in the large-r limit we ignore a small radius- 
dependent shift of lu. Then we find 



H 



vac+mat 



lu + u\ r fl 
-1 



(20) 



where 



Xr = A,- ^ = V2G F [n e (r) - n s (r)} ^ (21) 

encodes the net matter effect. Note that the 1/r 2 vari- 
ation here has nothing to do with the matter density 
profile, the effect of the latter appears through the n e (r) 
term. 



Henceforth we drop the explicit subscript r to denote 
the r-dependence of all quantities. Moreover, we drop the 
limits of integration which are always as above. In the 
small-amplitude case we have s UiU = 1. This simplifies 
in particular the diagonal Hamiltonian term which is 

K. u = 2 " + | / du 1 (u + v!) I dJ g u , tU , . (28) 

In the neutrino-neutrino term, the integral which involves 
J du'u' ... is a constant that does not depend on lu or u 
and therefore amounts to a shift of all frequencies, i.e. 
yet another rotating frame. Once more we can drop this 
term and arc left with 



^uj^u — 



lu + u (A + e/i) 



Here, 



€ = dudu) g^^ 



(29) 



(30) 



quantifies the "asymmetry" or "total lepton number" of 
the neutrino spectrum, normalized to the total v e flux. 
The EoMs arc then explicitly 



\d r Su, 



= [u + u(X + e/j,)] S u ,u 

- fj, / du' duj' (u + u 1 ) g u > u i S u 



(31) 



This is the linearized form of the EoMs and provides the 
starting point for the stability analysis. 

At this point it may be useful to recapitulate the el- 
ements that have gone into this analysis. Besides the 
small-amplitude approximation |S W , U | <C 1, we have 
taken the neutrinos to be far away from the neutrino 
sphere, R/r -C 1. At the same time, they have not 
yet reached the MSW resonance region, so that the or- 
dinary matter effect is large and the effective mixing an- 
gle in matter is small. We have also assumed that the 
vacuum mixing angle is so small that we may approx- 
imate cos 9 = 1, but it is trivial to accommodate an- 
other choice. We have assumed that the fast-rotating off- 
diagonal component of the Hamiltonian matrix caused by 
the mismatch between the mass and flavor directions av- 
erages to zero, so the only off-diagonal contribution of the 
Hamiltonian is provided by the neutrinos themselves. In 
numerical simulations, this fast-rotating component pro- 
vides the initial disturbance to kick-start exponentially 
growing modes. Here, however, we do not ask how the 
instability gets started, we only ask for the existence of 
exponentially growing modes. 



G. Eigenvalue equation 

The stability analysis determines if the small quanti- 
ties Su^u grow exponentially with r. This is achieved by 
writing S u u as 



-IQr 



(32) 



where both Q u ^ u and Q are in general complex numbers. 
A purely real solution for Q would imply that all modes 
precess with a common frequency. A complex solution 
£1 = 7 + in, with k > 0, would indicate an exponen- 
tially increasing S UtU , i.e., an instability. On the other 
hand, k < would indicate that the solution decreases 
exponentially toward the asymptotic solution S UlU = 0. 
In terms of Q u ,u, the EoM becomes 

(uj + uX - Q,)Qu,u = M / du' duj' (u + v!) g u > u i Q u ',u' , 



(33) 

where A = A + e/i. This may be looked upon as an eigen- 
value equation for Qui,u, which is a vector in the func- 
tion space on the uj-u plane, with the eigenvalue £7. The 
cigenfunction Q UlU implicitly carries a label £7 because 
usually for every £7 there exists a different Q^.u- Note 
that if Quj^u and its corresponding £7 satisfies Eq. (|3"3"|) . so 



does Q* u with the eigenvalue £7*. This implies that for 
each complex solution for £7 = j+in, there exists another 
solution £7 = 7 — in. Thus, the exponentially increasing 
and decreasing solutions always appear in pairs. 



III. SINGLE-ANGLE STABILITY ANALYSIS 
A. The consistency conditions 

Many important features of collective neutrino oscilla- 
tion phenomena can be understood in a simplified model 
where only a single angular mode u = uq is occupied, i.e., 
all neutrinos are assumed to be emitted from the neutrino 
sphere at a fixed angle relative to the radial direction. We 
therefore first perform the stability analysis with this as- 
sumption and later extend it to the multi-angle case. In 
the simplest schematic SN model, the neutrino sphere at 
radius R is pictured as a blackbody source without limb 
darkening, in which case the angular distribution is such 
that u is uniformly distributed on the interval < u < 1 . 
It is natural to represent this case in the single-angle ap- 
proximation by uq = 1/2. For the time being, however, 
we keep uq as a free parameter. 

Since u — uq, the term uX in the EoM corresponds to 
a common precession for all modes. We therefore can go 
to a basis rotating with frequency uqX, in which Eq. pip 
becomes 



idrSu, = u S u - 2u fi / dJ g u > S { 



(34) 



The single-angle approximation is then equivalent to say- 
ing that all the neutrinos feel the same refractive effect 
due to the other neutrinos. Requiring the solution to be 



of the form SV 



Qu 



-ifir 



gives 



(uj -ti,)Q u = 2«oM / dco' g^'Q 



(35) 



This is the single-angle form of the eigenvalue equation 
in Eq. (|55]l. 

For the l.h.s. of the eigenvalue equation to be indepen- 
dent of uj like the r.h.s., we must have 



Qlj oc 



1 



uj — n 



and therefore 



/i = 2uq / duj 



-n 



(36) 



(37) 



For an instability, this equation should have a complex 
root ft = 7 + ik. Then, splitting the equation into real 
and imaginary part, one obtains the two equations 



(2u fi) 1 = dug, 



duig u 



uj — 7 



(uj — 7) 2 + K 2 



(uj — 7) 2 + K 2 



(38) 

(39) 



Eqs. ([38)1 and (|39j) are the conditions that must be simul- 
taneously satisfied by 7 and the positive quantity k 2 . If a 
solution exists, we automatically have a pair of solutions 
fi = 7 + ±i|«|, of which one corresponds to an instability 
that will grow at the rate e' K '*. 

When an instability occurs, \Quj\ 2 takes the form of a 
Lorentzian centered at ui — 7 and half-width character- 
ized by k. Thus, the solutions for 7 and k tell us the 
range of w-modes which are significantly affected. We 
therefore present our results in the form of the parame- 
ters 7 and ft. Note that significant flavor transformations 
take place only if g u is significant in this range. 

In the form of Eqs. ([38]) and (|39|) . these results were 
previously derived in Rcf. fl9j where the full solution 
was provided, not only the small-amplitude expansion. 
The exponentially growing and shrinking solutions corre- 
spond to the "flavor pendulum" near its inverted position 
close to the beginning or end of a full swing. A purely 
real fl (no exponenti al g rowth) corresponds to the pure 
precession mode [HOI, HI- while Ec l s - <ES and ® 
are not new, the linearized analysis illustrates the origin 
of the Lorentz denominator in the eigenfunction Q u . 

A pure precession mode is described by a real eigen- 
value ft = 7. In this case Qu, becomes singular for lu = 7 
and the small-amplitude expansion is not self-consistent 
if g("f) 7^ 0. The full non-singular self-consistency re- 
lation for the pure precession mode, without the small- 
amplitude approximation, was provided in Ref. [15| and 
of course agrees in the appropriate limit. The linearized 
equations are useful to study the presence of instabilities, 
but not necessarily to study the pure precession solutions. 



B. A single spectral zero crossing 



The conditions in Eqs. (|3"8")l and (|3T))) allow us to un- 
derstand some of the stability features analytically. For 
example, for the integral in Eq. (|39[) to vanish, the inte- 
grand has to be positive in some parts and negative in 
the other. Thus, an instability requires the spectrum g^ 
to have a zero-crossing. The existence of such an insta- 
bility also requires the spectrum to cross from negative 
to positive values fl9| . 

During the accretion phase of SN evolution, powerful 
v e and v e fluxes are emitted, with a much weaker flux of 
other flavors. If we model this situation by a pure Fermi- 
Dirac spectrum of e-flavored neutrinos in the inverted 
hierarchy, the spectrum is positive for v e and negative 
for v ei providing for a single-crossed spectrum (Fig. [TJ. 
Therefore, instabilities are a generic feature of the neu- 
trino flux streaming from a SN core. 

Note that in the limit k — > 0, the integrand in Eq. (|3"9")l 
becomes g^irSiu — 7). Therefore, the consistency condi- 
tion implies 3(7) = 0. Thus in this limit, the instability 
starts near the zero crossing. Also, the r.h.s. of Eq. (|38|) 
gets a large contribution from the modes lu ~ 7, implying 
a small fj,. Therefore in the single-angle approximation, 
at small neutrino density one generically gets a narrow 




-1 1 

Frequency oj 

FIG. 1: Fermi-Dirac spectrum for v e and v £ , with Am 2 = 1, 
and T = 1 for both species. A small degeneracy parameter 
?7 = 0.28 provides an excess flux of v e over u e , which corre- 
sponds toe = 0.666. 



instability centered at a positive crossing. 



C. Normal vs. inverted hierarchy 

Recall that all our results have been obtained using in- 
verted hierarchy. Going to normal hierarchy corresponds 
to changing Eq. (|M|) to 



id r S u 



-w Su, - 2u [i / duj' g u > § ( 



(40) 



In terms of the solution S u of Eq. (|3~4")> . the solution of 
this equation is given by 



S u (n,g u ) = S*(n,-g u ) = S*(-n,g u ) 



(41) 



Since S and S* should have the same stability behavior, 
this implies that the stability conditions for normal hi- 
erarchy are the same as those for the inverted hierarchy 
with a change in the sign of g u . 

Formally we may implement this change of sign by 
keeping the same g^ and /i — > —fj,. In this sense one can 
show the solutions k(/j) and j(fi) on the same plot for 
both hierarchies by extending it to negative values of /i. 



D. Multiple spectral crossings 

A supernova emits neutrinos of all flavors, although 
with different fluxes and spectra. As g u is the difference 
spectrum between v e and u x , usually there will be ad- 
ditional spectral crossings other than the one at lu = 0. 
In the limit /i — > and concomitant n — > 0, there will 
be a solution Q^ centered on every positive crossing as 
already stressed in Ref. [19|]- For larger /j, there can be 
fewer solutions, but if there are several, they can coex- 
ist even when the different eigenfunctions are not well 
separated. 

In the single-angle case, our problem is equivalent to 
the reduced pairing Hamiltonian which is at the core of 



the BCS theory of superconductivity [32|. In this con- 
text, the same stability problem was recently investigated 
and it was shown how general analytic solutions can be 
constructed [331 ] . In other words, explicit large-amplitude 
solutions were constructed that correspond to several co- 
existing solutions with different H-values. The single-J7 
large-amplitude solution of Ref. [19[ is the simplest case 
of the class of multiple "normal soliton solutions." 



IV. BOX-SPECTRA EXAMPLES 
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In order to understand the behavior of the stability 
criteria it is useful to study explicit examples of single 
and multi-crossed spectra that can be solved analytically. 
The simplest approach is to represent a spectrum such 
as Fig. Q] in the single-energy approximation by one 5 
function for one average u (neutrinos) and one for neg- 
ative u> (antineutrinos) , or more such spikes to represent 
a multi-crossed spectrum. One then finds simple polyno- 
mial equations for tt with coefficients that depend on the 
frequencies and heights of the spikes. 

We here go one step further and represent the spectra 
by "boxes" of unit height, i.e. g w only takes the values 
0, ±1. This makes the integral in Eq. ([33]) once more 
analytically calculable and f2 becomes the root of a poly- 
nomial whose degree depends on the number of boxes. As 
will be seen, the simplified box spectra already bring out 
some important features of the stability of the realistic 
SN spectra. 

Scaling the height of all the boxes by a factor a corre- 
sponds to scaling the spectrum as g^ — > ag^ . The results 
for this may be obtained by simply scaling fj. — > a/i. For 
the sake of simplicity, in the analytic single-angle argu- 
ments we use Mo = 1/2. The results for any uq may be 
obtained by the scaling /i — > 2uq/i. The numerical results 
are given for various sample uq values. 



A. Two boxes 

Our first example is a schematic representation of the 
Fermi-Dirac spectrum of Fig. [T] We define our two-box 
spectrum as 



gu> = 



-a < w < 
< u < b 



(42) 



as shown in Fig. [51 We assume that all frequencies are 
normalized to some common frequency scale. The con- 
sistency condition in Eq. (j3"5)) yields 



n- 



(Q + a)(Q-b) 



V 



(43) 



where ij = e 1 ' A \ This corresponds to the quadratic 
equation 
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FIG. 2: Two-box spectrum with a = 1 and b — 1.6. Lower 
panels: 7(/i) and k(u), where broken curves are single angle 
with uo = 1, 1/2 and 1/4 (left to right). Thin solid (red) line: 
Multi-angle with uniform distribution < u < 1. Solutions 
for 7(u) are only shown where k(u) 7^ 0. 



from which the stability can be analyzed analytically. 
Note that < r/ < 1, while the limits /i->oo and /i->0 
correspond to rj = 1 and rj = 0, respectively. 
The roots of Eq. (gU) are 



-(b - a)r] ± a/(6 - a)V - Aabrjjl ^~q) 

2(1 - *?) ' l ' 

We have two complex conjugate solutions when the ar- 
gument of the square-root is negative, i.e. for 



< 77 < 77. 



Aab 



sync 



{b + af 



(46) 



(1 - i-/)Q 2 + (brj - ar])Q. + abi] = 



(44) 



Here Ji sy nc is the "synchronization strength," i.e. for a 
larger interaction strength the system is stuck in a stable 



position. If < 77 < 77 sync , wc find 



Three boxes 



7 = 



± 



(b-a)r} 

2(1-77) ' 

^Aab-q(l - rf) - (b - a) 2 r) 2 



2(1 -V) 



(47) 



These solutions arc shown in Fig. [2] for a spectrum with 
a = 1, b = 1.6, and different values of the single emission 
angle Uq. 

For r\ > rj sync the solutions for $7 are two different real 
roots so that k = 0. At the synchronization strength 
?7sync, w e find 



7sync 



lab 
Jb^a) 



(48) 



In our figures we only show the solutions with nonvan- 
ishing «, since when k = 0, the complex solutions 7 ± in 
do not exist and there are no instabilities. Therefore, the 
7(/i) curves are only shown where «(/i) 7^ 0. 

At the vanishing interaction strength /j, = 0, we have 
7 = k = 0. Had we put a spectral gap between the two 
boxes, the lower-/z cutoff point would be at a nonzero 
interaction strength fi > 0. 

For the completely antisymmetric two-box spectrum 
(a = b), Eq. (jll|) reduces to 



n 1 



abr) 



(49) 



which always has purely imaginary solutions. Thus, there 
is no synchronized behavior for such a spectrum, the fla- 
vor conversions take place at arbitrarily high values of \x. 
We will not further consider this special case that requires 
a fine-tuned spectrum. Moreover, a SN core always pro- 
duces an excess flux of v e over v e due to deleptonization. 
The spectrum shown in Fig. [2] has a positive zero- 
crossing. The case of a negative zero-crossing may be 
studied by multiplying the spectrum g u by a factor of 
— 1, i.e. g u — > —gu,. The consistency conditions then stay 
the same, with the change [i, — > — fj,, or 77 — > I/77. Equa- 
tion (14~4"1) then becomes 



(l-r))Q 2 -(b-a)Q-ab = 



(50) 



which always has real roots. Thus, the spectrum is stable 
as anticipated from our earlier arguments. 

To summarize, for the two-box spectrum in the single- 
angle approximation, the spectrum with a positive zero- 
crossing is always stable in the normal hierarchy In the 
inverted hierarchy (i) the spectrum with a negative zero- 
crossing is always stable, (ii) a completely antisymmet- 
ric spectrum with a positive zero-crossing is unstable at 
any value of /_*, (iii) an asymmetric spectrum with a pos- 
itive zero-crossing is stable above a threshold value of 
the interaction strength and unstable if the interaction 
strength is smaller. 



The natural next case is a spectrum with two crossings, 
of which only one can be positive, and so we expect at 
most one unstable solution for the eigenvalue equation. 
However, it was recognized earlier that such a spectrum 
shows additional features in that k(/j) need not cut off at 
large y, [281 ]. 

We represent this case by three adjacent equal- height 
boxes of the form (Fig. [3]) 



gu, = 



-a < lj < 
< uj < b 
b <uj < c 



(51) 



The consistency condition in Eq. (|35[) leads to the cubic 
equation 



(l-?7)£7 3 - (c - 2br] + arj)n 2 
-(b 2 - 2ab)rin - ab 2 ?] 







(52) 



It can have three real roots, or a single real root and a pair 
of complex conjugate roots. The latter case corresponds 
to instability. 

In order to study the stability at large /i, we look at 
the limit 77 = 1. The cubic term in Eq. (|52")l drops out, 
leaving us with the quadratic equation 



(a-2b + c)fl 2 + (b 2 - 2ab)Q + ab 2 = 

When the total lepton number is not zero (26- 
the solutions are 



n = 



2ab -b 2 ± yJb{Aab + b 2 - 4ac) 
2(a-26 + c) 



(53) 
•c^0) 

(54) 



If Aab + b 2 — 4ac > then both roots are real, implying 
that the system is stable. Wc then get the same cutoff 
behavior as in the two-box case, similar to the one shown 
in Fig. dJ 

On the other hand, if Aab + b 2 — Aac < 0, and writing 
Q = 7 ± in, we find 



7 : 



2ab - b 2 
2(a-2b + c) ' 



yJb{Aac -b 2 - Aab) 
2(a~2b + c) 



(55) 



Thus k reaches a non-zero asymptotic saturation value 
that is independent of /1. Likewise, 7 approaches a (i- 
independent asymptotic value, implying that the eigen- 
function Q u remains roughly centered on the nonvanish- 
ing spectral range. The system is thus unstable even at 
large [i values. This saturation behavior, shown in Fig. El 
was absent in the two-box scenario. This would allow fla- 
vor conversions to start deep inside the SN core. 

The spectrum in Fig. [3] has its central box positive. 
The corresponding spectrum with the central box neg- 
ative can be treated in the same manner, with the re- 
placement g^ — > — g u , which corresponds to /i — > — /j, or 
7/ — > I/77. In the 77 = 1 limit (large neutrino density) the 
behavior is then identical with the one described above. 
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unstable at arbitrarily large /i. For the inverted spectrum 
9u — > —Qu, the condition for stability is reversed. 

The three-box spectra thus show a cutoff or saturation 
behavior, depending on the details of the spectrum. The 
saturation behavior cannot occur for two boxes and is 
a new feature. The three-box spectra thus open up the 
possibility that for certain combinations of spectra and 
hierarchy, the neutrino ensemble is never stable. We will 
see, however, that this behavior is qualitatively modified 
by multi-angle effects. 
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FIG. 3: Three-box spectrum with a — 1,6 = 0.5, c = 0.6. 
Lower panels: y(fl) and k(/x), where broken curves are single 
angle with uo = 1, 1/2 and 1/4 (left to right). Thin solid 
(red) line: Multi-angle with uniform distribution < u < 1. 
Solutions for 7(^1) are only shown where k(/u) 7^ 0. 



When the total lepton number vanishes, Eq. ([5 
duces to 

(1 - ?y)0 3 - (26 - a) (1 - 7])n 2 - (b 2 - 2ab)7]fl - ab 2 T] = . 

(56) 
One can examine if it has three roots of the form £7i , Vl? = 
7+ik, O3 = 7 — in in the limit rj = 1. Clearly in this limit, 
17 1 = ab/(2a — b) is a real finite root. Then the sum of 
the roots is fli + 2j = 26 — a, a finite number, implyting 
that 7 is finite. On the other hand, the product of the 
roots, 17i(7 2 — k 2 ) = —ab 2 r\j(\ — rj) is a large negative 
number. This is possible only if k is a large number, and 
f2i > 0, i.e. 2a — b > 0. Thus, the stability in this case 
is determined by the sign of (2a — b). If this is negative, 
then there is a cutoff behavior, i.e. the system is stable 
in the 77 = 1 limit. If this is positive, then the system is 



C. Four boxes 

A four-box spectrum comes closest to representing a re- 
alistic SN spectrum which typically has three zero cross- 
ings as explained earlier. Representing it by four adjacent 
boxes, we define the spectrum as 



-1 


—a < u < — b 


+1 


-b < uj < 


-1 


< w < c 


+1 


c < LU < d 



9lj = 



The self consistency condition is 

(n + b) 2 (n~c) 2 



(57) 



(n + a)n 2 (n-d) 

and can be written as 



(58) 



(1 - f])^ + (2b -2c- 077 + dr])n 3 

+ (b 2 - Abe + c 2 + adri)fl 2 

+(2bc 2 - 2b 2 c)VL 



7.2 2 

be 







(59) 



This quartic equation can have zero, one or two pairs of 
complex solutions, so this is the first explicit case with 
the possibility of two simultaneous unstable solutions. 

We are interested in the high-density behavior where 
77 = 1. In this limit, the quartic equation reduces to 
a cubic equation with real coefficients. Hence we are 
guaranteed at least one real root and we can have at most 
one pair of complex conjugate solutions, i.e. at most one 
solution for k 2 . Its existence is determined by the value 
of the discriminant 

A = & 2 c 2 [-276 2 c 2 (a-2& + 2c-rf) 2 
+ 32b(b-cfc(-a + 2b-2c + d) 
- 366(6 - c)c(-a + 2b-2c + d) (b 2 - Abe + c 2 + ad) 



A(b-c) 2 (b 2 -Abc + c 2 + adf 
4(6 2 -46c + c 2 + ad) 3 ] . 



(60) 



If A > 0, all the roots are real and the system is stable. 
A typical scenario is shown in Fig. 01 At low [i values, 
there are two instabilities, one of which shows a cutoff be- 
havior, i.e. it vanishes for fi greater than a certain value. 
The other instability, with 7 « 0, survives for arbitrarily 
large values of ^t, with k showing a saturation behavior as 
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FIG. 4: Four-box spectrum with a — 1.5, b = 1.0, c = 0.6 and 
d — 1.6. Lower panels: 7(^1) and k(aO, where the solution 
with larger 7 is termed number 1. Broken curves are single 
angle with uo = 1, 1/2 and 1/4 (left to right). Thin solid 
(red) line: Multi-angle with uniform distribution < u < 1. 
Solutions for 7(/z) are only shown where n([i) ^ 0. 



in the 3-box case. Under the transformation g u — > —g u , 
fi — > — /i, t] — > 1/77, the instability condition at 77 = 1 
does not change, so that at large neutrino densities the 
saturation solution exists in either hierarchy. 

The four-box spectra can give rise to two instabilities 
at the same time, or only one, or none at all, depending 
on the spectral details and strength of /i. At most one 
of them can show saturation behavior at large density, 
which then exists for both hierarchies. 



V. SPECTRA WITH TAILS 

The ultimate goal of our investigation in the context of 
SN physics is to understand in which regions of the star 
the neutrino stream may show self-induced instabilities. 
In agreement with the previous literature we have found 
that in some cases, and notably in the generic two-box 
example, there is a value /J, sync such that for pi > /i sync 
the system is stable. Upon closer scrutiny, however, one 
finds that the stable regime docs not exist for a realistic 
spectrum. The main difference between the Fermi-Dirac 
spectrum of Fig. [1] and the two-box spectrum of Fig. [2] 
is that the former has nonvanishing tails for large |w|, 
corresponding to the infrared part of the spectrum. In 
other words, g^ ^ everywhere except at the spectral 
crossing point. Given a positive spectral crossing, the 
only solutions with k — > can be those centered on the 
spectral crossing and corresponding to fj, —> 0. Purely 
real solutions for il then do not exist for large fi and the 
system is always unstable. 

One can easily illustrate this point by a numerical solu- 
tion of the eigenvalue equation for the Fermi-Dirac spec- 
trum of Fig. [I] We show 7(/i) and k(h) in Fig.[5l assum- 
ing single angle with uq = 1/2. The result is qualitatively 
similar to the two-box example, except that at large /z the 
curve k([i) has a tail. 

In a SN core, the system as a function of radius sweeps 
from very large ^-values (of order 10 6 in our normaliza- 
tion) to zero, but the system is nowhere stable. Still, the 
usual onset of bipolar oscillations happens approximately 
at the /j,- value around the would-be synchronization value 
where the "hump" and the "tail" of the k(^) curve join, 
in Fig. [5] around [i ~ 10-11, as confirmed in numerous 
numerical SN examples. 

The explanation is two-fold. At large \i the system is 
always unstable, but k(//) is very small and the exponen- 
tial growth will not get anywhere on the available time 
scale (or rather, radial distance scale), in particular as 
the center of the resonance, 7(/i), is constantly shifting 
as a function of the decreasing fi. Second, the center of 
the resonance is far away in the spectral tail (the infrared 
energy modes), so only a narrow range of infrared modes 
is affected. It deserves mention that "small k" simultane- 
ously means two things: The exponential growth is slow 
and the resonance width is narrow. 

Evidently, then, a stability analysis alone is not enough 
to understand the consequences in a realistic SN. One 
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FIG. 5: Eigenvalue solution for the Fermi-Dirac spectrum of 
Fig. [l] shown by the brown (solid) line. As a dotted (black) 
line we show the same result if the Fermi-Dirac spectrum is 
cut off for |w| > 4. Bottom panel: (p). 



(P> = 



Loo dujp u 



(62) 



We can speak of an effective cutoff behavior if (p) as a 
function of \i quickly drops to very small values in a nar- 
row range of jx, while being appreciably large at smaller p. 
In the bottom panel of Fig. [5] we show (p) for the Fermi- 
Dirac example. Indeed, for /x ~ 11, this measure drops 
by orders of magnitude and we would have an effective 
cutoff — above this transition range, only a minimal range 
of infrared modes experiences self-induced conversion. 

Another approach is to artificially cut off the infrared 
part of the spectrum. The curve 7(/J.) in Fig. [5] has a 
distinct kink exactly in the region where the the hump 
and tail of the k(/i) curve join. If we cut off the Fermi- 
Dirac spectrum and set it to for \u>\ > 4, the k(/j,) 
curve follows the original one almost exactly, except that 
the tail is cut off at /U syn c ~ 11 and the system becomes 
perfectly stable for larger \x values (Fig. [5]). Therefore, 
the simple picture of a cutoff behavior is actually very 
good except for a minor infrared correction and our study 
of box spectra indeed provides a useful picture of the 
instability behavior. 

The saturation solutions that can arise for spectra with 
two or more crossings (three or more boxes) are very dif- 
ferent. The saturation solutions have a k that is compa- 
rable to the width of the spectrum and remain centered 
on the main part of the spectrum. Therefore, such an 
instability cannot be removed by truncating the infrared 
part of the spectrum. These instabilities are only tamed 
by multi-angle effects. 



VI. MULTI- ANGLE STABILITY ANALYSIS 



A. The consistency conditions 



needs a more dynamical approach that takes account of 
the available time scale as p, sweeps from large to small 
values, and one needs to consider the range of modes that 
are actually affected by the resonance. 

One rough way to approach the latter point is to con- 
sider the maximum total flavor conversion that can be 
achieved by a solution $7 on a given spectrum, or equiva- 
lently, the average conversion probability over the entire 
spectrum when the flavor pendulum has moved from the 
inverted to the normal position. From the exact large- 
amplitude instantaneous solution [19( one finds that the 
maximum conversion probability of a mode of frequency 
lo is 



(61) 



(bJ — "f) 2 + K 2 



The linearized stability analysis in the single-angle 
case shed new light on the self-consistency conditions as 
originating from a simple eigenvalue equation, although 
the conditions themselves had been found previously to- 
gether with the full solution independently of the small- 
amplitude expansion. However, the single-angle approxi- 
mation is not justified in the SN context except that nu- 
merically single-angle simulations and multi-angle simu- 
lation sometimes, but not always, yield similar results. 
Therefore, based on the eigenvalue equation ()33|) . we 
now extend the linearized analysis to the multi-angle 
case. Besides the small-amplitude expansion, this equa- 
tion also uses the large-distance approximation where it 
is assumed that the angular divergence of the neutrino 
radiation is small. 

The crucial step is to realize that the r.h.s. of Eq. ([55)1 
is of the form A + Bu where A and B are expressions 
that do not depend on cither uj or u. So we are led to 
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the ansatz for the form of the eigenfunction 

a + bu 

ui + uX — \l 



(63) 



where a and b are complex numbers. Inserting this ansatz 
provides 



a + bu = n I du duj' g u 



(u + u')(a + bu') 
u' + u'X - 



(64) 



To understand better the structure of this equation, we 
define the integrals 



/„ = / dudujg UtU 



uj + uX — Q 
Then our eigenvalue equation becomes 



a + bu = (x 



(ah + bh) + (al + bh)u 



(65) 



(66) 



If this is supposed to be true for every u we need to match 
the coefficients of the linear u polynomial on both sides 
separately. We can then write this in matrix form 



/' 



h h 
h h 



(67) 



This has the form of an eigenvalue equation for a 2 x 2 
matrix. This equation has nontrivial solutions if 



det 



h h 
h h 



/' 







or explicitly 



/'• 



1 = h ± y/hl- 



2 ■ 



(68) 



(69) 



This is the multi-angle counterpart of our single-angle 
eigenvalue equation of Eq. (|3"T)) . 

We introduce once more f2 = 7 + in, and to split this 
equation into its real and imaginary parts we write the 
integral expressions in the form 



where 



•hi 



I71 "n 1 l-"-n j 

u + uX — 7 



(70) 



duj dug u . u u n 



(uj + uX — 7) 2 + k 2 



K n = / dujdug U:U u r ' 



[bJ + llX — 7) 2 + K 2 



(71) 



Inserting this into Eq. (|69| and equating the real and 
imaginary parts separately yields the two real equations 



(Ji - A'" 1 ) 2 = K 2 + J0J2 - K K 2 , 
_!, J K 2 + K„J 2 



(J x - /i- 1 ) 



2A'i 



(72) 



These two equations can now be taken as our consistency 
conditions, analogous to Eqs. (f3"51) and (pl9"|) in the single- 
angle case. 



Note that for a given angular mode u = uq, the quan- 



tity \Q U 



is a Lorentzian, centered at ui = 7 — uqX. 



Thus, the range of w-modes that are affected by the insta- 
bility is different for different angular modes, as opposed 
to the single-angle case. 

In order to solve for 7 and k satisfying the consistency 
conditions, one may consider \i and A as independent 
quantities, keeping in mind the constraint A = A — e/i. 
This allows us to find one real equation from which fi 
is eliminated. As in the single-angle case, this equa- 
tion provides a /i-indcpendent relation between 7 and 
K and thus the set of all possible eigenvalues. The set 
of all points in the fl = 7 + in plane fulfilling this re- 
quirement represents the "root locus diagram" of our 
self-consistency relations. The other equation is of the 
form n~ l = (integral expressions) and allows one to cal- 
culate for any allowed Q = 7 + in the corresponding /i by 
quadratures alone. 

The energy and angular distributions of SN neutrinos 
are not independent of each other. On the other hand, 
in schematic models one may assume that the angular 
distribution is independent of energy and independent 
of flavor. In particular, the approximation of a com- 
mon blackbody neutrino sphere for all flavors implies 
9ui,u = 9u, independently of u. In this case the u in- 
tegration in the expressions for I n can be performed ex- 
plicitly, considerably simplifying the numerical inversion 
of the self-consistency equations. 



B. Matter vs. neutrino background 

Multi-angle effects on the stability of the spectrum 
may be interpreted in terms of the separate effects of 
the two terms on the r.h.s. of the EoM in Eq. (jST]) . The 
second term represents the effects of Pantalconc's off- 
diagonal refractive index caused by neutrino-neutrino in- 
teractions. It allows for the instability in the first place 
and in addition can lead to self-induced multi-angle deco- 
herence [26|,[27|]. However, this decoherence effect arises 
after some pendular oscillations and would not be visible 
in the small-amplitude expansion. Decoherence arises in 
a periodic system performing many revolutions and dif- 
ferent modes dephasing relative to each other. The un- 
stable solution in the small-amplitude expansion is not 
periodic, but exponentially growing, and thus we can 
not study decoherence effects. However, we can study 
the multi-angle modification of the instability which is 
caused by A in the first term on the r.h.s. of Eq. (|3"i"T) . 

If Pantalconc's flavor off-diagonal refractive effect 
would not exist, i.e. the last term on the r.h.s. in 
Eq. (|3ip would be missing, then the first term on the 
r.h.s. specifics that every mode S u ,u freely processes 
around the weak-interaction direction with frequency 
w e g = uj + (X + efi)u. We recall that 



A = \/2GF{n e — Tie) 



2r 2 



(73) 
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as given in Eq. (|2Tj) . represents the net matter effect in 
the co-rotating frame. Note that we have kept only the 
matter term causing an angle-dependent spread of w e ff, 
while the main matter effect, causing a common preces- 
sion with frequency \/2GF(n e — rig), has been removed 
by going to a rotating frame. 

The flavor-diagonal refractive effect caused by the neu- 
trino background, in our equations represented by en, 
plays a perfectly analogous role to ordinary matter. In- 
deed, it is the quantity A = A + e/i which appears in the 
EoMs, and hence in the multi-angle analysis. Though 
we define the neutrino background contribution in terms 
of the radial flux densities and not their number densi- 
ties, at large distances these two quantities are identical. 
Therefore, we may write 



B 2 



en 



\/2G F [(?v -np e )- {n Vx - n 5 J] r-o 



2r 2 



(74) 



The matter background term A and the neutrino back- 
ground term e\x are completely symmetric, if one takes 
into account that the x contribution is missing among 
charged leptons simply because of the absence of charged 
x leptons since both /i and r leptons are too heavy to be 
present in the SN context. 



C. Multi-angle suppression of the instability 

A large matter effect in the form of a large value of 
A = A + e/i can suppress the self-induced instability in 
numerical SN simulations. This effect was predicted in 
Ref. |3(| and has been repeatedly tested, most recently in 
the numerical studies of Ref. [24[ . The reason invoked for 
this suppression was the large dispersion of w e a = u>-\-Xu 
that is caused if u is spread over the unit interval and A 
is much larger than the direct spectral spread of w. It 
was then argued that this multi-angle matter suppression 
of the instability would be noticeable if the matter term 
were much larger than the neutrino background term, 
i.e. if A ^ ii. The linearized analysis allows us to show 
analytically that when A ^ [i, the collective oscillations 
are completely suppressed. 

Let w max be the largest value of |w| where g^^ is signif- 
icant. When A ^> u max , the quantities /„ in Eq. ([53)1 are 
suppressed due to the factor of Xu in the denominator. 
The cancellation of large values of Xu by il is possible 
only in a narrow range of u values, which will give only 
a small contribution to the integral. Equation (|69p can 
then be satisfied only at fi ~ A. Therefore, the system 
should be stable for A S> \x, i.e. when the net effects 
of ordinary matter dominate over those of the neutrino- 
neutrino interactions. 

A seemingly different multi-angle suppression affects 
the "saturation mode" of collective oscillations that can 
occur in spectra with at least two zero crossings for very 
large u fSecs. lTVBl and lrV C[) . This mode was first discov- 
ered in an analytic toy model |28j and later re-discovered 



in the context of a more realistic SN example [29( . More- 
over, both papers argued and tested numerically that 
the spread of uj c s caused by the neutrino-neutrino term 
would suppress the saturation mode. This suppression 
effect, however, is just another aspect of the impact of 
A in the EoMs. The matter and neutrino backgrounds 
together cause a common multi-angle effect by dispers- 
ing w G ff. This may be seen as follows. 

In the single-angle approach and for spectra with e not 
very small, the consistency relations imply that k cannot 
be much larger than a typical spread Aw of the frequency 
spectrum. Indeed, for k ^> Aw, Eq. (|38| would not be 
satisfied for any finite asymmetry e. Assuming that k 
and thus the resonance width remains of order Acj in 
the multi-angle case, it is clear that a large dispersion 
of w e s shifts most modes away from the resonance. In 
this way, multi-angle effects would indeed suppress self- 
induced flavor conversions at large A. We will see that 
the presence of a large A in the EoMs will never increase 
k significantly, although it can suppress it considerably. 
In this sense the reasons given here for the multi-angle 
matter suppression indeed apply. 

For certain spectra, notably the generic two-box case, 
the single- angle approximation predicts a cutoff /^ S ync- 
The instability cannot form at large values of \x since 
the consistency conditions in Eqs. (|38|) and (|39|) cannot 
be satisfied for any complex Q, and the system is stable. 
For fi < u sync , however, pendular oscillations commence. 
The multi-angle analysis does not change this behavior 
qualitatively. The value of ^ sync in the multi-angle case 
is some average of the ^t S ync values for different uq values, 
as can be seen in Fig. [5] 



D. Normal vs. inverted mass hierarchy 

In the single-angle approximation, we saw that as far 
as the stability analysis is concerned, analysis of normal 
hierarchy is the same as that of the inverted one, except 
for a change of the sign of u. This is true also in the multi- 
angle scenario, except that one also needs to change the 
sign of A. Indeed, normal hierarchy changes Eq. (|3"Tj) to 

id r Su, u = [— w + u(X + en)] S u<u 

- n / du duj' (u + u) g u > u ' S u >, u ' ■ (75) 



The solution of this equation can be given in terms of the 
solution S u u of Eq. (J3TJ) as 



bu>,u{n, *,9u 






(76) 



Since S and S* should have the same stability behavior, 
this implies that the stability conditions for normal hi- 
erarchy are the same as those for the inverted hierarchy 
with a change in the sign of g u or \x (not both at the 
same time), and an additional change in the sign of A. 
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VII. MULTI-ANGLE EXAMPLES 

Wc now study the stability conditions including multi- 
angle effects for the box spectra considered earlier as well 
as for a realistic SN example. It turns out that the multi- 
angle effects modify the single-angle results in significant 
ways. We always take the emission to be uniform over 
< u < 1. The integrals J n and K n can then be ana- 
lytically calculated for the box spectra. The expressions, 
however, are not very illuminating, and we do not give 
them here. 



A. Two boxes 

The two-box spectrum and the single-angle eigenvalues 
j(ii) and k(li) for different choices of the emission angle 
Mo were shown in Fig. [5J In addition, we show in Fig. [2] 
the numerical result for 7(/x) and k(ii) in the multi-angle 
case, assuming the absence of matter so that A = e/z and 
all effects are caused by neutrino-neutrino interactions 
alone. As far as the stability behavior is concerned, the 
multi-angle effects may be interpreted as some (compli- 
cated) average of the single-angle effects with different 
uq values. In particular, the value of k at any value of 
\x lies within the range of k for different uq values in the 
single-angle approximation. In other words, multi-angle 
effects leave the system qualitatively unchanged. 

For sufficiently dense matter the picture changes con- 
siderably. In Fig. [6] we repeat the single-angle results 
for uq = 1/2 (thin black dot-dashed curves) and the 
no- matter multi-angle results (red solid curve). We also 
show curves for A = 5 and 30. Qualitatively the 7(/x) and 
k(li) curves are shifted to larger values of it, the spectra 
being stable below some threshold value and above some 
cutoff value. The interval of unstable /x values remains 
roughly the same, except that it is shifted to larger li val- 
ues. The shift is approximately such that the unstable 
domain arises for it, ~ A as predicted earlier. 

The presence of matter alone does not stabilize the sys- 
tem, it shifts the instability domain to larger \x values. 
Indeed, a nonzero A can allow for an instability even at 
the values of ii where e/x alone could not have generated 
an instability. In the li-X plane, the system is unstable 
along a strip roughly following A ~ /i, and stable outside. 
In a SN we sweep from large to small ii values as a func- 
tion of radius, and at the same time from large to small 
A values. If everywhere along the path A 3> li, we would 
be outside the instability strip. 

However, even if this is not the case, the system could 
still be stabilized. We also have the effect of dispersing 
the spectrum relative to the resonance condition. In a 
spectrum of given spectral width Acj, the resonance de- 
nominator u> + uX — il will be on resonance only for a 
narrow range of u values if A is large. Therefore, even if 
the system is not stable, most modes will be off-resonance 
in analogy to our discussion of spectral tails. 




FIG. 6: Multi-angle eigenvalues 'j(li) an d k(aO for the two- 
box spectrum in Fig. [5] Black (dash-dotted) thin line: single- 
angle (uo = 1/2). All other lines: multi-angle. Red (solid): 
no matter (A = 0). Brown (dashed): A = 5. Purple (dotted): 
A = 30. The 7(/z) curves are only shown where k(/x) 7^ 0. 



A detailed understanding of the multi-angle stabiliza- 
tion of the neutrino flux streaming from a SN core as 
studied most recently in Ref. [24j requires more detailed 
scrutiny because the described effects of shift and disper- 
sion both affect a realistic SN neutrino flux. 



B. Three boxes 

We now turn to the three-box example of Fig. [3] which 
was chosen such that we observe a "saturation mode:" 
The system is unstable with an asymptotic value of k 
for arbitrarily large li and the center of the resonance 
7 remains centrally located in the spectral range. Here 
the multi-angle matter effect causes a dramatic modifi- 
cation of the single-angle results because the curve «(//) 
is completely suppressed above a critical /x-value, i.e. the 
spectrum is actually stabilized. 

In Fig. [7]we repeat the single-angle case with uq = 1/2 
and show the multi-angle case with vanishing matter so 
that A = tfi. The previous saturation effect is completely 
suppressed. As we increase A, the instability domain is 
shifted to larger li values roughly linearly with A as ex- 
plained earlier. In other words, the system now behaves 
similar to the two box case. 
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FIG. 7: Multi-angle eigenvalues 7(a0 and k(m) f° r the three- 
box spectrum in Fig. [3] Black (dash-dotted) thin line: single- 
angle (uo = 1/2). All other lines: multi-angle. Red (solid): 
no matter (A = 0). Brown (dashed): A = 5. The "/(fi) curves 
are only shown where k(aO 7^ 0. 



C. Four boxes 

The main novelty of the four-box spectrum is that it 
can display two simultaneous instabilities. One of them 
can be of the saturation type, similar to the three-box 
example, details depending on the exact choice of pa- 
rameters. Multi-angle effects for A = and thus A = e/x 
are shown in Fig. 01 The results are entirely in line with 
the expectations from the previous examples in that the 
saturation mode gets suppressed at large \x. We shall not 
repeat the analysis with the four-box spectra here. How- 
ever the realistic SN spectrum in the next section can be 
seen to be close to a four-box spectrum, and is observed 
to display the same expected behavior. 



D. A realistic SN spectrum 

We now apply our stability analysis to a realistic SN 
spectrum that is motivated by Duan and Friedland's 
recent numerical study of multi-angle suppression j29l |. 
They used v e , u e , and v x spectra with assumed Fermi- 
Dirac form. The temperatures were taken to be 2.1, 3.5, 
and 4.4 MeV, respectively, and the degeneracy parame- 
ters 3.9, 2.3 and 2.1. The average energies are then 9.4, 
13.0 and 15.8 MeV. The ratios of number fluxes are taken 
to be 1.3 : 1.0 : 1.5, where the total v e flux is normalized 



to unity as per our convention introduced in Sec. Ill Al 
We use Am 2 = (50meV) 2 in order to convert the energy 
scale to uj, which we show in the units of km -1 . 

The spectrum is shown in the top panel of Fig. [5] It 
has three spectral crossings and thus compares with a 
four-box spectrum, although the deviation from zero of 
the left-most part is almost invisible. (For ui < —0.8, the 
value of g u is small and negative.) The main feature of 
interest here and in Duan and Friedland's study is the 
presence of a "saturation mode" in which the instability 
in the single-angle treatment exists for arbitrarily large 
/j. Such a feature requires at least two crossings, so in 
this regard the visual impression of this being a "three 
box" example is actually the main point. 

In single angle (thin dot-dashed black line) we find a 
dominant "saturation mode" as expected (corresponding 
to 72, K2) and in addition a "cutoff mode" (correspond- 
ing to 71, «i) which always has very small values for k(/x). 
The existence of this mode is a residue of the not per- 
fect disappearance of g^ at w ;$ —0.8, i.e. of the small 
"fourth box." If we were to truncate the spectrum at 
the left-most crossing, this mode would disappear. At 
very small \x, these two instabilities are centered near the 
positive spectral crossings as expected. Even the "cutoff 
mode" does not completely disappear for large fi due to 
the spectral tails. 

Next we include multi-angle effects without matter 
(A = 0), roughly corresponding to the study of Ref. [29[ 
where the matter density was small compared to the neu- 
trino density. The red solid curve represents the case with 
A = efj,. We see that at large \i, corresponding to small 
radii in the SN, the saturation mode is indeed suppressed. 
It does not vanish completely, owing to the presence of 
spectral tails, in contrast to the earlier box examples. 
The center of the saturation mode, 72 — ue/i, continues 
to be in the central part of the spectrum for u s=s 1. The 
center of the cut-off mode, 71 — ue^j,, on the other hand, 
is actually pushed into the tail of the spectrum. 

Next we add a significant amount of matter, repre- 
sented by A = 5 (brown dashed curves). There now a 
lower ^-threshold appears for both, the cutoff mode as 
well as the saturation mode. This is a consequence of the 
complete suppression of collective oscillations at A ^$> (x, 
predicted in Sec. I VI Cl Thus, our analytical prediction is 
vindicated in this scenario of a realistic SN spectrum. 

The main result of Ref. j29(, obtained in the regime 
A ~ fj,, was that multi-angle effects stabilize the spectrum 
in deep layers of the SN where in single-angle strong bipo- 
lar oscillations would have occurred. Our stability anal- 
ysis provides two reasons for this behavior. The value of 
At(/i) is strongly reduced, although the spectrum is always 
unstable as behooves a spectrum with tails. Moreover, 
the dispersion of effective oscillation frequencies due to 
the presence of uX in the resonance denominator moves 
most modes off resonance. Given ^t(r) and A(r), one can 
predict the exact SN radius where serious bipolar oscil- 
lations would kick in for the multi-angle case, with the 
stability analysis. 



16 






-2 


-1 




Frequency a> 


1 


2 


60 










, - 


4(1 










■ 


20 

f- 
-20 
-40 






72 




: 












-60 










""-•- = 



50 



100 150 



200 




200 




100 
A' 



150 



200 



FIG. 8: Realistic SN spectrum as described in the text and 
identical with Ref. [2Jj]. Lower panels: j(fJ-) and «(/i) for 
the two instabilities, where 1 denotes the one with smaller 7. 
Mode 2 is the saturation mode. Black (dash-dotted) thin line: 
single-angle (uo = 1/2), otherwise multi angle. Red (solid): 
no matter (A = 0). Brown (dashed): A — 5. 



VIII. MULTI-ANGLE INSTABILITY IN 
NORMAL HIERARCHY 

A perfectly symmetric spectrum shows self-induced 
kincmatical dccoherence among angular modes in ei- 
ther hierarchy, or in the language of our present study, 
whether the zero-crossing of g^ is positive or nega- 
tive [23. The same is found for spectra with a small 
asymmetry e< 1 [27j . Kinematical decoherence requires 
the spectrum to be unstable in the first place, so these 
findings imply that in the multi-angle case, an instabil- 
ity does not require a positive zero-crossing of g u . This 
surprising result is easily verified with our method. 

Let us assume that A = 0, as would be the case for 
an antisymmetric spectrum (e = 0) and in the absence 
of matter. We further make the simplifying assumption 
of a universal angular distribution, i.e. g UjM = g u f u , with 
J duf u = 1 . Let us introduce the notation 



G= duj 



uj — Q 
Then in the multi-angle case, we have 

I n = G f duu n = G(u n ) . 
Jo 

The eigenvalue equation is then 

M- 1 = h ± VV2 = \(u)±(u 2 ) 1 / 2 



(77) 



(78) 



(79) 



The stability analysis then corresponds to determining 
the solutions £1 = 7 ± in for the equation 



M 



K±G 



(80) 



where K± = (u)±(u 2 ) 1/2 - 

In the single-angle approximation, (u) — u , (u 2 ) = Uq, 
so that fC+ = 2uq and /C_ = 0. The latter equation has 
no solution, since it would require \jT x = identically. 
One may then write 



^single 



2u () G 



(81) 



This may be solved to obtain the values of 7 and k in 
the single-angle approximation. Clearly, /x^ n le has the 
same sign as G, so that G is positive for a spectrum with 
positive crossing. 

We can now write the two values of /i -1 in multi-angle 
analysis that correspond to the same values of 7 and k: 



/'+ 



= ICjlG 



II- 



1 = tC-G 



(82) 



Thus, there are two values of /j, in multi-angle analysis 
that correspond to the same (7, n) in the single-angle 
analysis. The multi-angle stability behavior at these two 
(i values can then be understood in terms of the single- 
angle stability behavior at ^ s i ng i e . 

Note that /C_|_ > and /C_ < 0. Therefore, /i+ is 
positive, and it corresponds to just a scaled value of 
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Msingie with the same (inverted) hierarchy. The single- 
angle study is indeed a good proxy for the multi-angle 
case in such a scenario. On the other hand, /i_ is nega- 
tive, and corresponds to the normal hierarchy. For spec- 
tra with (u 2 ) 1 / 2 ~ (u) one finds that \n~\ S> \(J>+\, so 
the "wrong hierarchy" case shows the same instability at 
much larger densities. 

The exact mapping between the single- and multi-angle 
cases is not possible for the more realistic case of A being 
small but nonzero, but of course we expect a qualitatively 
similar behavior. 

In Ref. [231 ^ was found that kinematical decoherence 
is suppressed for sufficiently large e, which in normal hier- 
archy means that the spectrum is stable and the "positive 
crossing" rule for an instability is back in force. These 
studies did not include the presence of matter, so A = e/i 
with a sufficiently large e implied that A was not very 
small. In this sense, the stability in normal hierarchy for 
sufficiently large e looks like yet another case of multi- 
angle suppression of an instability. 



IX. CONCLUSIONS AND OUTLOOK 

The refractive effect of neutrino-neutrino interactions 
can cause a flavor instability among the different energy 
and angular modes, exchanging the flavor content while 
leaving unchanged the overall flavor content of the entire 
ensemble. In the context of SN neutrinos, the focus so far 
has been on impacts of collective flavor oscillations on the 
neutrino signal of the next nearby SN or for the r-process 
nucleosynthesis. Many analyses have been performed un- 
der the simplifying assumption that all neutrinos feel the 
same refractive effect due to the other neutrinos (single- 
angle assumption). However it has also been found that 
multi-angle effects can strongly modify the answer and 
can typically suppress flavor conversions that would oc- 
cur in the single-angle treatment. In this sense, a more 
pressing question is to understand the stability condi- 
tions for a dense neutrino stream, and to understand if 
the usual neglect of flavor oscillations in those dense SN 
regions is justified where a Boltzmann treatment of neu- 
trino transport is necessary. 

We have studied the conditions under which instabil- 
ities can occur in the neutrino stream in the two-flavor 
framework. While the equations of motion are nonlin- 
ear, the run-away solutions can be found by a standard 
mode analysis of the linearized system, i.e. in the small- 
amplitude expansion. We have shown how the stability 
question is reduced to a straightforward eigenvalue equa- 
tion in the co-rotating frame. The existence of complex 
eigenvalues directly corresponds to an instability. In the 
multi-angle case, besides the small-amplitude expansion, 
we have linearized the equations in the small angular di- 
vergence of the neutrino flux that is applicable at a large 
distance from the source. In the single-angle case, our 
conditions agree with the previous literature, whereas in 
the multi-angle case, our results are new. 



In the single-angle case and using schematic energy 
spectra (box spectra in the u> variable), the eigenvalue 
equation can be solved analytically. We identify two 
generic instability cases. The "cutoff mode" exists for 
a finite range of neutrino densities, or equivalently, a 
finite interval of effective neutrino-neutrino interaction 
strengths /Uo < M < Mi- I 11 particular, the system is 
stable above a cutoff value fi\, also termed the "syn- 
chronization" or "sleeping top" regime. This behavior 
is generic to a single-crossed (or two-box) spectrum and 
represents the traditional flavor pendulum. On the other 
hand, there is the "saturation mode" where the insta- 
bility approaches an asymptotic value for large \x. This 
mode requires at least two crossings (three boxes) and 
had been previously found in the literature. When it ex- 
ists, the saturation mode exists for both hierarchies. For 
multi-crossed spectra, several modes can co-exist, for ex- 
ample in the four-box case a saturation mode and a cutoff 
mode or two cutoff modes. 

An additional case occurs for neutrino spectra with a 
vanishing asymmetry, i.e. with no net lepton number flux. 
The run-away rate n grows indefinitely as a function of /z. 
This actually is a special case of the cutoff mode (cutoff 
at infinite /x) and does not seem to be of practical interest 
in the SN context. 

Realistic neutrino spectra are not described by boxes 
and have tails as a function of cj = \Am 2 \/2E, corre- 
sponding to infrared energies. Such spectra do not have a 
cutoff and a synchronization regime does not exist. How- 
ever, for large /z the instability moves entirely into the 
infrared regime and has no practical impact. There re- 
mains an effective cutoff density where the amount of 
flavor exchange between modes drops by orders of mag- 
nitude within a narrow range of neutrino densities. Still, 
the existence of a synchronization regime is only an ap- 
proximate concept. 

Multi-angle effects modify the instability in several 
ways. For example, while the single-angle approximation 
predicts an instability only if there is a positive zero- 
crossing of the spectrum g^, the multi-angle analysis im- 
plies that spectra with a small asymmetry (a small net 
lepton number flux) and in the absence of matter are un- 
stable in both hierarchies, whatever be the sign of the 
zero-crossing. The solutions for both hierarchies can be 
mapped onto each other. 

The main multi-angle modification of the instability is 
caused by a background of matter, where here neutrinos 
themselves also contribute to the matter effect. The sat- 
uration mode is "fragile" and, for box spectra, develops 
a cutoff. (For realistic spectra with tails, there is no ex- 
act cutoff, however, as mentioned earlier.) The neutrino 
background alone is enough to achieve this effect. 

For both the saturation and cutoff modes, a large mat- 
ter effect shifts the domain of instability zzo < M < A f i to 
larger values of (i. The amount of shift is approximately 
equal to the matter potential A, apparently without re- 
ducing the typical k values within the instability domain. 
An additional consequence is that the affected range of 
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modes shrinks for increasing A. Even if the run-away 
scale k is not reduced, only a narrow range of angular 
modes is on resonance if A is large. For realistic spectra 
with tails, k(h) does not vanish for large fi, i.e. there is 
not an exact cutoff. On the other hand, there is a lower 
threshold /j.q ~ A below which the spectrum is completely 
stable. In other words, multi-angle effects introduce a 
range < \l < [Aq ~ A where the spectrum is stable even 
when such a stable range does not exist in the single- 
angle approximation. This behavior corresponds to the 
multi-angle suppression of instability at A 3> /U, which we 
motivate analytically through our linearized treatment. 

One may further investigate these findings in the con- 
text of realistic SN situations where neutrinos stream 
from regions of large density to vacuum, sweeping 
through \x values from very large to zero. The central fre- 
quency 7 of the instability then sweeps through different 
spectral ranges and the growth rate k depends on both 
the neutrino and matter density. Effective flavor conver- 
sion can be suppressed in different ways: The growth rate 
k can vanish when /i falls into the multi-angle-implicd 
stability domain, or it can be too small on the available 
distance scale, or the instability can be located mostly 
in the infrared and affect only an ignorablc part of the 
spectrum. 

We have considered neutrinos streaming from a source 
and have made the approximation of small angular di- 
vergence that is appropriate at a large distance from the 
source. Our approach can be easily adapted to a homo- 
geneous ensemble with a nontrivial angular distribution 
that evolves in time rather than as a function of distance 
from a source. In this case no expansion in the angular 
variable is needed. 

We have only considered trivial angular distributions 
where all neutrino flavors are assumed to be emitted from 
a common blackbody sphere without limb darkening, cor- 



responding to a box spectrum in our u variable. One can 
extend our analysis to more complicated angular distri- 
butions that could involve, for example, spectral cross- 
ings in the angular variable. As advocated by Sawyer, 
such crossings cause a "multi-angle instability," but its 
practical relevance remains to be investigated. 

We have made the usual assumption of azimuthal sym- 
metry of the neutrino flux around the radial direction 
away from the source. In the long-distance limit of small 
angular divergence of the neutrino field, one can of course 
expand the EoMs in terms of small variables that de- 
pend on two angles. It remains to be explored if cylin- 
drical symmetry enforces unphysical solutions, just as the 
single-angle approximation sometimes enforces unphysi- 
cal solutions, and if there are new instabilities or impor- 
tant modifications that could be relevant in the SN con- 
text. Based on schematic models using a small number of 
modes it has been suggested that violations of cylindrical 
symmetry could have a major impact |l3j |. 

A rich phenomenology of neutrino flavor conversions is 
bound to emerge from further investigations into the col- 
lective effects due to neutrino-neutrino interactions and 
ordinary matter background, in the complete multi-angle 
framework. 
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